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ABSTRACT 

Low-ionization-state far-infrared (FIR) emission lines may be useful diagnostics of 
star-formation activity in young galaxies, and at high redshift may be detectable from 
the ground. In practice, however, very little is known concerning how strong such line 
emission might be in the early Universe. We attempted to detect the 158 pm [C ii] 
line from a lensed galaxy at z = 4.926 using the Caltech Submillimeter Observatory. 
This source is an ordinary galaxy, in the sense that it shows high but not extreme star 
formation, but lensing makes it visible. Our analysis includes a careful consideration of 
the calibrations and weighting of the individual scans. We find only modest improve¬ 
ment over the simpler reduction methods, however, and the final spectrum remains 
dominated by systematic baseline ripple effects. We obtain a 95 per cent confidence 
upper limit of 33mJy for a 200 kins -1 full width at half maximum line, corresponding 
to an unlensed luminosity of 1 x 10 9 L 0 for a standard cosmology. Combining this 
with a marginal detection of the continuum emission using the James Clerk Maxwell 
Telescope, we derive an upper limit of 0.4 per cent for the ratio of T[cii]/-^fir in this 
object. 
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1 INTRODUCTION 

Low-ionization-state far-infrared (FIR) emission lines play 
an important role in cooling star-forming regions, and 
they allow us to infer the flux of rest-frame ultravio¬ 
let (UV) photons. Star-forming galaxies have large dust 
masses which obscure UV radiation, so we rely on the 
IR and longer wavelengths to probe these star-forming re¬ 
gions. Studies have shown that [Cn], the ground-state fine- 
structure line of singly ionized carbon ( 2 Pz /2 — > ~P 1 / 2 , A = 
157.7409 pm), is the domi nant co oling line in gas-rich star- 
forming regions llDalgarno fe M cCravlll972l : ICrawford et alJ 
ll985l : ISt,acev et allE)fuli ~ It is thus expected to be a probe 
of star formation in young galaxies. In some nearby galaxies 
the [C 11 ] line accounts for a s much as 1% of the far-IR lu- 
minosity dStacev et al .1 [ToOll : iNikola et al • 1998 : Mochizukil 
2000). On the other hand, in some more distant ultralumi- 
nous infrared galaxies (ULI RGs) results from ISO indicated 
much lower frac tions <|Malhotr a et al.ll200li IContursi et alJ 

l2002l : lLuhman et alfeOO.L At 2 ~ 5 we know very little con¬ 

cerning the physical conditions within star-forming galaxies. 
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Metallicity, UV flux, dust content and geometry could all 
be radically different from that for more nearby (and hence 
cosmically older) systems. In the absence of much empirical 
data, even upper limits on [C 11 ] luminosities are therefore 
worthwhile. 

[ C 11 ] is the brightest line in our Galaxy dMizutani et all 
ll994Tl . but its wavelength makes it impossible to observe 
from the ground. At high redshift, where the line is red- 
shifted into the submillimetre, there is great promise for us¬ 
ing this line to probe early galaxy evolution, or even to fin d 
the first objects (e.g. ISta.rkjll99A ISuginohara et, al III 99Tb 
Until instruments with much wider back-ends and greatly 
improved sensitivity become available, progress can only be 
made slowly, by targetting individually promising objects. 

The z = 4.926 lensed galaxy, detected as a red arc in 
the c luster CL 1358+62 (referred to as ‘Gl’ in lFranx et alJ 
Il99?l hereafter F97), magnified by a factor of 5 — 11, rep¬ 
resents a particularly promising target for observing [C 11 ] 
from the ground. Its redshift places the line at an accessible 
frequency (~ 320 GHz), though not ideal since it is on the 
edge of an atmospheric water line. The galaxy has a reason¬ 
ably high inferred star formation rate (~ 36 Mg yr -1 , F97), 
and unlike quasars observed at similar redshifts, the [C 11 ] 
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line is expected to be close to the systemic redshift of the 
galaxy. Moreover, although rest-frame UV colours indicate 
significant reddening, there is little submilli metre cont i nuum 
detec ted, with a 3<r upper limit of 4 mjy (Ivan der Werf et alJ 
2001, but see Section [T|. 

In the following paper we discuss the analysis and in¬ 
terpretation of our data. Section |21 briefly describes the ob¬ 
serving programme. Section [3 describes the specifics of our 
analysis procedures. The results are presented in Section 0 
and discussed in Section^] 


2 OBSERVATIONS 

Object G1 was targetted with the Caltech Submillimeter 
Observatory (CSO), located at 13 300 feet on the summit of 
Mauna Kea in Hawaii, over three nights in 1998 January. 
The telescope is a 10.4 m //0.4 submillimetre antenna with 
an alt-az mount. 

The atmospheric transmission during the observing run, 
as measured by the CSO r meter at 225 GHz, was very good. 
The optical depth ranged from T225 — 0.03 to 0.07, corre¬ 
sponding to optical depth at 320 GHz of 0.32 to 0.55. 

2.1 Frequency Centre 

F97 determined the redshift of G1 from the detection of 
Lya emission at 7204 A, giving z = 4.926. They note, how¬ 
ever, that the Ly a line is assymetric and that the measured 
Sill absorption line is blue shifted by ~ 300kms _1 . These 
features can be explained by an outflow model, whereby an 
expanding shell of neutral and i onized materi al abs orbs the 
blue side of the Ly a line (e.g. iLeaueux et a DE995;). We as¬ 
sume that the [C 11 ] emission would be somewhere between 
the Lya and Sill lines, and therefore centred the detector 
shifted Av « 200kms _1 from the Lya centre, or z = 4.922. 
This is a reasonable assumption, as this shift is typical of 
linewidths observed in local galaxies. The [C 11 ] emission line 
is redshifted to uq = 320.939 GHz for the expected centre of 
our source. 


2.2 Instrumentation 

Based on the calculations of the previous section, it is clear 
that the CSO 345-GHz receiver is the appropriate detector. 
This receiver is a single side-band SIS mixer with 1 GHz IF, 
and acousto-optical spectrometer (AOS) backends. For our 
observations the signals from the mixer were sent to three 
spectrometers simultaneously, with 50-, 500-, and 1500-MHz 
bandwidths. 

From the F97 detection of Sin in CL 1358+62-G1, we 
assume the [C 11 ] line is ~ 200kms _1 wide, corresponding 
to 200 MHz at our observing frequency. We therefore con¬ 
centrate our efforts on data from the 1500-MHz receiver, use 
the 500-MHz data as a check and ignore the 50-MHz data. 

On the first night the receivers were centred on the nom¬ 
inal line centre, 320.939 GHz. In an attempt to avoid confus¬ 
ing any systematic detector response with the shape of the 
emission line, the receiver was shifted by Av = 129.8 km s _1 
for the second and third nights. 

The secondary mirror chop throw was set to 60 arcsec at 


Table 1. Details of observations. Row 1 shows the line centre 
of the receiver on each night. The rest of the table displays the 
number of usable scans taken on each night in each detector. 


Night: 


1 

2 

3 

Receiver centre (GHz): 

320.939 

320.800 

320.800 

1500 MHz 

Source scans: 

142 

236 

42 


Calibr. scans: 

42 

71 

10 

500 MHz 

Source scans: 

142 

236 

42 


Calibr. scans: 

41 

75 

10 


1.123 Hz on the first night. It was apparent that after obser¬ 
vations of one night we had reached near the instrumental 
noise limit, so we set more conservative chop parameters, 
40 arcsec at 0.7 Hz, on the second and third nights in an 
attempt to probe systematic effects. By changing the chop 
throw, we additionally decrease the probability of inadver¬ 
tantly chopping on to a nearby source. 

2.3 Pointing and Calibration 

The telescope was pointed at a — 13 h 59 m 39 f !0, S = 
62° 30'47”. This position was determined by visually com¬ 
paring the Hubble Space Telescope (HST) image presented 
in F97 wit h co ordin ates of cluster members published by 
iLunpino et all (1991), as F97 did not publish coordinates 
for Gl. This process is accurate to within a few arcsec, 
which is well within the CSO beamsize at this frequency 
(~ 20 arcsec). Pointing calibrations were performed occa¬ 
sionally throughout each night, using the 12 CO (3 —> 2) line 
of IRC 10216. 

Temperature and frequency calibrations were also per¬ 
formed several times throughout each observing session. 
Hereafter, ‘calibration scan’ refers to the temperature cali¬ 
brations, i.e. spectra of a hot load, while ‘source scan’ refers 
to spectra of our source, Gl. 

I 11 all, ~ 450 usable source scans were acquired, simul¬ 
taneously in each of the 500 and 1500 MHz detectors. Exact 
numbers are given in Table Q The 1 o RMS per raw scan is 
~ 40 — 80 mK. 


3 ANALYSIS 

CSO spectroscopic data are usually analysed using the 
CLASS (Continuum Line Analysis Single-dish Software, 
iBuisson et aljfl995 ) software package. In the expectation 
that any detection would be at very low signal-to-noise ra¬ 
tio, we developed our own analysis procedures in order to 
maximise our control over the reduction process. In this sec¬ 
tion we describe our methods of calibration and reduction. 

3.1 Calibration 

CSO data are calibrated in real time by referring each source 
scan to the previous hot load calibration scan. Visual in¬ 
spection of the 500- and 1500-MHz calibration scans show 
apparently unstable baselines, so we investigate whether a 
more careful calibration will result in lower noise. 

We compare the mean value of each calibration scan to 
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Night One Two Three 




Figure 1. Trends in data as a function of time: (a) is the mean 
value of each calibration scan for both the 1500- (squares) and 
500-MHz (triangles) spectrometers, measured in K; (b) is the vari¬ 
ance of the source scans for both spectrometers, in K 2 ; (c) is the 
optical depth measured at 225 GHz with the CSO r meter. Note 
the correlation between all three panels. The vertical dashed lines 
indicate the divisions between the three nights. 

the variance of the following source scans and to the atmo¬ 
spheric opacity, T 225 - We find that the 500- and 1500-MHz 
calibration means and the source scan variances all correlate 
well, and that they track r (see Fig. 0. 

We propose three possible methods of calibrating the 
data: (i) for each source scan, use the most recent calibra¬ 
tion scan (the CSO default); (ii) naively average together all 
calibration scans within each night; and (iii) interpolate each 
calibration channel in time across each night. The data are 
thus reduced, as described in Section m using each of the 
three calibration methods. Taking the variance in the final 
reduced scan as a measure, we find that method (ii) is clearly 
inferior to (i), while (i) and (iii) are not significantly differ¬ 
ent. We choose method (iii) as the most elegant solution. If 
baseline variations did not dominate the error budget, we 
might expect method (iii) to yield the best results, and we 
would advocate its use in general for low signal-to-noise ratio 
data. 


3.2 Data Reduction 

The data are co-added within each night. Firstly, a zeroth- 
order baseline, fit to the wings of the spectrum, is subtracted 
from each source scan in order to minimize the effects of 
baseline drift. The data are then combined as a weighted av¬ 
erage in each channel. Because of small shifts in the velocity 
centres and spacing throughout each night, we interpolate 
the data to a common set of velocity bins before averaging. 

The summed scans for each night are binned into 
roughly 40 kms -1 wide channels, and the nights are aver¬ 
aged, weighted by the variance in each bin. We rebin the 


Figure 2. The co-added source spectrum, using data from the 
1500-MHz detector. Error bars indicate statistical errors, calcu¬ 
lated as the standard deviation of the data in each bin. The spec¬ 
trum is dominated by baseline ripples. The solid curve is the 
95 per cent upper limit for the [C Ii] emission line centred at 
v = Okms -1 (320.939GHz) and with FWHM= 200kms -1 . 

data to six different bin sizes and centres, and find no sig¬ 
nificant difference in the final result due to the choice in 
binning. 

Finally, the co-added spectrum is converted from an¬ 
tenna temperature to flux density. The beam size and effi¬ 
ciency for the 345- GHz receive r are 24.6 arcsec and 74.6 per 
cent, respectively, llKoo liH. The resulting spectrum is 
shown in Fig. [^] 


4 RESULTS 

It is clear from Fig. [5] that the error bars are much smaller 
than the structure in the final spectrum. The systematics in 
the instrument prevent integrating down to such low noise 
levels. We re-estimate the error bars by finding the best-fit 
Gaussian and scaling the errors so that = N, where N is 
the number of degrees of freedom in the fit. 

A full four-parameter Gaussian fit will not converge to 
a meaningful result, so we instead fix the line centre and 
width to the a priori expected values of vo = Okms -1 and 
full width at half maximum (FWHM) = 200kms -1 (as de¬ 
scribed in Section 12.11 . and calculate the baseline as the 
weighted average of the data outside ±2cr of the line cen¬ 
tre. We find a 95 per cent upper limit of 33mJy using the 
1500-MHz receiver; the 500-MHz receiver gives consistent, 
although less constraining results. To convert from specific 
flux S v to luminosity L, we first integrate over frequency, 

5 = f S„dv = 1.06 ST ak x FWHM for a Gaussian of 
height SJ( eak and width FWHM = uo(Av/c), then multi¬ 
ply by 47where Dl is the luminosity distance for a 
given cosmology. For a standard cosmological model with 
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Hm = 0.3, Ha = 0.7, and Ho = 70kms -1 Mpc -1 , we find 
Dl = 45.8 Gpc and A[cn] = 5.4 x 10 9 L®. 

We also reredu c e exis iting SCUBA observations at 
850 pm llHolland et allflOOflli to find an estimate of the total 
FIR luminosity, Afir. An 850-u.m upper lim it of S„ < 4 mjy 
was reported bv lvan der Werf et al.l ( 2001 1 based on 4ks of 
pointed photometry data. We combine 15 ks of jiggle-map 
data (~ 2.5 mjy rms) with the phot ometry data (~ 1.5 mjy 
rms) using the methods described in lBorvs et al . (1200 30 and 
find a 2.1a peak 2arcsec from Gl, well within the 15-arcsec 
beam. Although not highly significant, and in a blank field 
survey we would not consider a 2a peak a detection, this is a 
pointed measurement of a known source, so we consider this 
a measurement rather than simply an upper limit. We can 
then calculate the ratio A[®n]/AFiR using Gaussian proba¬ 
bility distributions to describe both Ljcii] and Afir. Inte- 
grating a grey-body model with Td us t = 40 K and /3 = 1.5 
llCombesil999l : lKlaas et al.l200lll from 8 to 1000 pm, we then 
find that the SCUBA measurement Ssso = (2.8 ± 1.3) mjy 
corresponds to A fir = (2.4 ± 1.2) x 10 12 L®, assuming the 
same cosmology as above. Integrating over the joint prob¬ 
ability, we find a Bayesian 95 (99) per cent upper limit 
A[cii]/Afir < 0.4 (1.1) per cent. 


5 DISCUSSION 

Typical values of A[cii]/Afir are thought to be ~ 0.1 — 
1 per cent, but there is recent evidence to suggest that it 
can be lower, especially for ULIRGs, which are thought to 
make up a large fraction of the high-redshift submillimetre 
population. We find A[cii]/Afir < 0.4 per cent at 95 per 
cent confidence for our partic ular target. While this does 
not rule out the ratios found bv IStacev et al.l Il99 1). it does 
confine the ratio for this particular high-redshift galaxy to 
the lower end of the range. 

Based on their lensing model, F97 determine that the 
source Gl is magnified by a factor of 5 —11. Taking the lower 
end of this magnification range, we find an unlensed lumi¬ 
nosity of A[cn] < lx 10 9 L®. We can use the absolute Lpn] 
measurement to place an upper limit o n the star f orma tion 
rate (SFR) occurring in the source. iBoselli et all f 2002tl re¬ 
late SFR to A[cn] by assuming an initial mass function, 
with slope a = 2.35 between 0.1 and 100M®, and using 
their [C il] - Ho? relation determined from observations of a 
sample of Virgo galaxies, finding 

SFR = 1.73 x 10“ 6 x (Afcni/A©) 0 788 Meyr -1 . 

Assuming this relation holds for all galaxies at all red- 
shifts, and that the source is not differentially lensed, we 
find SFR < 21 M® yr -1 . However, we note that feoselli et all 
l2002t) claim the relation only holds up to 10 10 5 L® and that 
the SFR is uncertain by up to a factor of 10, so we should not 
place too much confidence in the estimate. As a comparison, 
F97 find SFR ~ 36M®yr _1 , based on near- i nfrared pho- 
tometry and the star formation model bv lBruzual fc Gharlotl 
( 1993h . Hence the SFR implied by our [Cn] upper limit is 
loosely consistent with that inferred from the optical. How¬ 
ever, given the approximations and assumptions made here, 
it is hard to draw any firm conclusions. 

iLuhman et alJ (2003) find a deficit of [Cn] (A[cn] < 


10~ 3 Afir) in a sample of ULIRGs, and claim a negative cor¬ 
relation between A[ C h]/Afir and A^(60 pm)/AM 1 00 pm ) 
for galaxies of all types (see Fig. 5 in llinhma.n et inEm). 
In order to compare our limit with sources in this plot, we 
calculate 

Af IR = 1.26 x 10 12 x [2.58 A„(60 H + A„(100 pm)] , 

with A, evaluated at rest-frame wavelengths (assuming 
a modifi ed blackbody) in WHz -1 and Af IR in W (see 
iHelou et alJIlOSsh . We find L[c n]/L' Fm < 4.4 x 10~ 3 and 
A„(60 pm)/A„(100 pm) = 0.88, which is consistent with the 
published relation for normal star-forming galaxies. Never¬ 
theless, firmer upper limits at similar levels for a sample of 
high redshift star-forming galaxies will help to determine if 
ther e is a dif ferenc e in the role of [C n] in the early Universe. 

iBolatto et alJ (’200j) recently reported an upper limit to 
A[cii]/Afir inaz = 6.42 quasar, observed with the JCMT. 
They find A[cn]/AFiR < 5 x 10~ 4 . While this limit is lower 
than that published here, we note that a quasar is not a 
normal galaxy, and this source is ~ 200 brighter than ours, 
with Abol ~ 10 14 L®. Our measurement, because of lensing 
amplification, has allowed us to place a useful limit in a 
single high-redshift star-forming galaxy. Clearly, this is a 
difficult programme, and more sensitive measurements of 
both the [C il] line emission and the FIR are required to 
significantly improve the constraint on A[c n i/AFiR. 
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